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METHOD FOR DETERMIimi S A VA U O E GIVEW ■ re -MKB i R R R m' P ARAMETERS 

A SYSTEM 



The present invention relates to a method for deter- 
mining the value to be given to a set of so-called specific 
parameters of a system based on the values of a set of so-called 
system measurement parameters. 
5 Such a method may be used to control various systems 

such as a character recognition system, an electric cotponent 
failure diagnosis system, a system for evaluating a transport 
cost.^ 

Fig. 1 shows an exanple of a system xising such a 
10 method. A car 1 attenpts to follow a truck 2, the car and the 
truck being models moving on a planar obstacle- free surface. The 
car is equipped with a camera 3 and with an autonomous control 
system. The control system has the function of defining at 
regular intervals, for exanple, every 100 ms, the direction and 
15 the speed that must be taken by car 1, as well as the inclina- 
tion of the camera so that the truck permanently remains in the 
camera's field of vision. 

At a given time to# the control system estimates where 
the truck will be 100 ms later, at a time t^. In the case where 
20 the truck effectively moves as expected, it will be at the esti- 
mated position Pe(ti), at the center of image 4 that the camera 
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will take at time ti- However, genersd-ly, the truck will be at a 
measured position Rm(ti) different from Pe(ti). Position Eta(ti) 
is shifted in the x and y direction with respect to Pe(ti) 
according to a horizontal shift Cx and a vertical shift Cy, 
5 Horizontal shift Cx indicates whether the truck has moved more 
to the left or to the right. Vertical shift Cy indicates whether 
the truck has accelerated or slowed down. 

Fig. 2A is a side view of car 1 moving at speed v 
provided with camera 3, and shows the vertical inclination angle 

10 of the camera's axis. Fig. 2B is a top view of car 1, and 

shows horizontal inclination axis of the camera's axis, as 
well as rotation angle cn^ot wheels of car 1. 

Fig. 3 shows a possible configuration of car 1 and of 
truck 2, at time t^, truck 2 being at position Pm(ti) of Fig. 1. 

15 The truck wo\ild thus seem to go rightwards (assuitption a) . The 
truck may also go forwards and accelerate (assiinption b) . The 
truck may aJ-so start back off to the left (assumption c) . It 
will be considered hereafter that the truck has but these three 
possibilities and is only likely to go to one of the three esti- 

20 mated positions Pe(t2)a, Pe(t2)b, and Pe(t2)c, 

The car control system must then decide whether the 
Ccu: must tcdce direction d^^ enabling joining truck Pe(t2)c or 
whether it must take direction d2 enabling joining the truck at 
position Pe(t2)a or Pe(t2)c. Since the three possibilities 

25 Pe(t2)a, Pe(t2)b, and Pe(t2)c are a priori equiprobable , the 
selection of direction d2 seems to be the most judicious since 
it covers a larger number of possibilities, and it will be 
considered hereafter that d2 has been selected. 

The control system must then choose to increase or 

30 decrease speed v of car 1. The image analysis enables thinking 
that the truck has accelerated since it is at the top of the 
image. A first decision could be to require an acceleration of 
the Ccu:. However, the control system has chosen to go in direc- 
tion d2 and it is possible for the truck to smoothly turn 

35 rightwards to reach to the position of truck Pe(t2)a. Knowing 
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this probability, it then seems judicious not to accelerate too 
much to avoid colliding with the truck. 

The control system similcurly defines angles ay and 
to be given to the camera according to the previously- taken 
5 decisions and to the estimate of the future displacement of the 
truck. 

The control system of car 1 must be able to define the 
new values of parameters v, a^ot' ^ ^h' with a small-size 
and low-cost calculator and memory. It is further necessary for 

10 the control system to take decisions quickly, every 100 ms. The 
control system must also process a large number of measurement 
parameters c^, Cy and of specific or, more specifically, of 
control parameters v, o^ot' ^ able to take efficient 

decisions that enable following the truck whatever its trajec- 

15 tory. 

The forming of a control system requires previously 
defining the interdependencies between parameters. Thus, in the 
above-described exanple, the selection of horizontal inclination 
angle ah and of rotation angle a^ot depends on the measured 

20 horizontal shift Cx; the selection of vertical inclination angle 
Oy. and the selection of speed v depend on the mecisured vertical 
shift Cy. Further, angle ah and angle a^t ^® dependent from 
each other, otherwise the truck risks coming out of the field of 
vision of the camera, whereby the car can then no longer follow 

25 the truck. Similarly, speed v and angle are dependent from 
each other, angle Oy having to be adapted to the speed and vice 
versa. 

A model of the joint probability distribution of all 
the system parameters is defined based on a set of independent 
30 probability distributions defined for each of the previously 
identified interdependencies. Thus, probability p(v, a^ot' 
ah# Cxf Cy) for a given combination of values of all the parame- 
ters to be possible and clever, can be defined for the above- 
described set by the following fonttula: 
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p(v, Orot' ^h' ^X' ^) = 

p (Cx) P (Cy ) p (tth/cx) P (Ov/cy ) p (ttrot/oth) P ( v/a^) ( 1 ) 
where p(cx) is the prdbability for horizontal shift Cx to have a 
given value, p(Cy) is the probability for vertical shift Cy to 
5 have a given value, p(ah/cx) is the probability to have a 
horizontal inclination angle knowing shift Cx, pCov/cy) is 
the prdbability to have a vertical inclination angle knowing 
shift Cy, p{a^ot/^h) probability to have a rotation angle 

<^rot knowing angle a^, and pCv/o^) is the probability to have a 

10 speed v knowing angle o^. 

Siirple nan- conditional prci>abilities, such as p(Cx)# 
may be defined by an analytic function, for exairple, 
p(Cx) = (absolute value of k/cx) where k is a nomnalization 
constant. Similarly, conditional probabilities, such as 

15 p(ah/cx) , niay be defined by a family of ancdytic functions, that 
may for exanple be gaussian functions centered on value Cx- 

Simple probabilities, p(Cx)/ or conditional probabili- 
ties, p(ah/cx)f can also be calculated from a database indexing 
the number of times when a horizontal shift Cx or a couple of 

20 values (tthf Cx) has been observed, for example, in a trai n i n g 
phase. 

The joint probability distribution model, the analytic 
functions, and the databases are memorized. After an image has 
been taken at a time t^, the control system must decide of the 

25 value to be given to each control parameter v, oc^^ot' ^h 
based on the values noted down for each measurement pairameter at 
time ti, Cxi, ^yi- system first selects a couple of 

values (cXj-of ' then, after, the other couple of values ((Xy, 
a^) . Only the selection of a couple (a^-ot' will be described, 

30 the selection of a couple (o^, a^) being performed identically. 
Probability pCorot'^/^xi'^yi) for the selection of a couple 
(a^-ot'^) to be pertinent knowing shift values Cxi and Cyi can be 
calculated according to the following formula: 

P(an)t.v/c,i,Cyi) = ^ ZP(^'^«)t'^vAh,C,i,Cyi) (2) 
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where Z is a normalization constant. 

There are two sorts of existing control systems capa- 
ble of selecting a couple (ttrot'^) based on the probability 
distribution of the couples. The first ones select a couple 
5 (ttrot''^) directly from formula (2) . The second ones construct, 
prior to the selection of a couple (o^ot''^^ ' ^ representation of 
the probability distribution of the couples. 

Among the first systems, so-called optimization 
methods which search the maximum probability and select the 

10 corresponding couple can be mentioned. Now, it may be preferable 
to select a couple which does not exhibit a probability maximum. 
For exanple, in the case where the probability distribution 
exhibits several maximums for two couples (o^rot, I'^^l^ 
(^rot,2'^2)' it ^ pertinent to select a couple having a 

15 speed intermediary between v^ and V2 and a rotation angle 
between OLj^ot,i and CLj^ot,2' Other methods such as the METROPOLIS 
method (see "Monte Carlo simulation and numerical integration" 
by J. Geweke, 1996), perform a random drawing of a couple 
directly from fonmila (2) . Such methods iitplement calculation 

20 algorithms, for the most part very complex and xise powerful 
calculators. Such methods further require a significant calcula- 
tion time, inconpatible with certain uses such as that described 
in Fig. 1. 

Among the second systems, certain methods use a so- 
25 called tabular representation of the probability distribution. 
This tabular representation consists of memorizing, for each 
couple (a^-ot'^) ' probability p(cXi;.ot''^/^xi'^yi) obtained for the 
shifts Cxi, Cyi nieasured at time ti. The selection of a couple 
may then consist of taking the maximnjm probability couple, or of 
3 0 performing a random drawing of a couple from the memorized data. 

In the case where the studied parameters are nxamerous, 
or when a parameter can take many values, the calculation time 
of the probabilities of each couple (a^ot/V) quickly becomes 
very long. The size of the memory used to store the tabular 
35 representation may also quickly become prohibitive. 
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A representation of the probability distribution, 
known as a "gaussian mixture", consists of modeling the distri- 
bution by a set of gaussian functions, each gaussian function 
being defined from a maximum probability value. Previously, an 
5 optimization method is \ised to identify the couples (Orot'^) 

maxiimim probabilities. The coi?)les are then distributed in 
groups containing more or less couples according to whether the 
couples have or not values a^ot ^ close to the identified 
maxiimams- A probability value is calculated for each group, the 

10 probability values forming gaussians around the maximum 
probabilities. The selection of a couple is then performed by 
random drawing or search of the couple of maximum probability. 

This method enables having a more or less accurate 
representation according to the available memory space. However, 

15 an increase in the accuracy of the representation requires a new 
distribution of the couples and a new calculation of the 
probabilities of each group. Further, the modeling of the 
probability distribution by a set of gaussian functions is not 
adequate for all systems. 

20 An object of the present invention is to provide a 

method for determining the values to be given to one or several 
specific parameters of a system knowing one or several meaisxire- 
ment parameters, in the case where the nimiber of parameters is 
very large and/or in the case where some of the parameters may 

25 take a great n\mber of values. 

Another object of the present invention is to provide 
a method for determining the values to be given to all the 
specific parameters of a system within a time interval that can 
be very short. 

30 Another object of the present invention is to provide 

such a determination method using a memory of variable size and 
possibly very small. 

Another object of the present invention is to provide 
such a determination method using a sinple calculation device. 



k 



To achieve these objects, the present invention 
provides a method for determining the value to be given to a set 
of specific parameters of a system based on the values of a set 
of measiirement parameters of this system, where each of the 
5 parameters can take a finite nxmiber of values, the system being 
associated with a means for providing a probability value for 
any combination of values of the specific parameters, said 
probability value being all the greater as the selection of the 
considered ccatibinatian is pertinent knowing the value of the 
10 measurement parameters, the method conprising the steps of: 

- noting down the value of each measurement parameter; 

- constructing a tree-shaped representation of the 
probability distribution of all the possible combinations of 
values of the specific peurameters corresponding to the noted 

15 down values, the set of conibinations, forming a first branch, 
being divided into several subsets of conibinations, forming 
second branches, each subset gathering conibinations having close 
specific parameter values, where each second branch can 
similcurly divide into several third branches and so on, a prob- 

20 ability value being assigned to each branch, this probability 
value being that obtained for one of the combinations of the 
considered branch or for one of the combinatians of one of the 
branches from which the considered branch originates; 

- selecting according to a predefined selection crite- 
25 rion one of the connbinations of values of the specific parame- 
ters based on the representation of the previously- constructed 
tree-shaped probability distribution. 

According to an inplementation mode of such a method, 
the branches resulting from the division of a same branch are at 
30 the number of two and contain the same number of combinations, 
the first branch dividing in two second branches, where each 
second branch can divide in two third branches and so on. 

According to an inplementation mode of such a method, 
the division of a branch in two branches conprises the steps of: 



a) selecting a ccnibination different frcm the combina- 
tions having already been xased to define the probability value 
of the existing branches and calculating the probability of this 
selected conibination; 
5 b) dividing the so-called "parent" branch containing 

the selected conibination in two so-called "child" branches; and 

in the case where the selected combination and the 
"parent" coinbination xised to define the probability value of the 
parent branch belong to the same child branch, assigning to the 
10 two child branches the probability value of the parent branch 
and dividing the child branch containing the selected combina- 
tion by resuming the method at step b) , this child branch 
becoming the parent branch, and 

in the case where the selected combination and the 
15 parent coinbination do not belong to the same child branch, 
assigning the probability value of the selected combination to 
the child branch containing the selected combination and assign- 
ing the probability value of the parent combination to the other 
child branch. 

20 According to an inplementation mode of such a method, 

the selection criterion consists of selecting one of the combi- 
nations exhibiting the maximum probability. 

According to an inplementation mode of such a method, 
the selection of a combination consists of iitplementing the 

25 recvirsive method ccitprising the steps of: 

a) randomly selecting a number p ranging between 0 and 

1; 

b) calculating the siam of the probability values 
assigned to the two so-called child branches resulting from the 

30 division of the first branch, and calculating for each child 
branch a new probability value equal to the ratio between the 
probability value assigned to this child branch and the calcu- 
lated sum; 

c) defining two contiguous probability intervals 
35 between 0 and 1, the first interval being associated with a 
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first child branch, the second inter\/eJ. being associated with 
the second child branch, the first interval ranging from 0 to 
and including the probability value of the first child branch 
and the second interval ranging from this probability value to 
5 1; 

d) identifying in which interval number p is to be 
found and selecting the child branch associated with this inter- 
val, and 

in the case where the selected child branch ramifies 
10 into other branches, resuming the recursive method at step a) , 
the first branch being replaced with the selected child branch, 
otherwise 

e) selecting one of the combinations of the selected 
child branch. 

15 According to an implementation mode of such a method, 

the selection criterion consists of selecting one of the combi- 
nations having a probability value which is predetermined or 
ranging between two given probability values. 

According to an inplementation mode of such a method, 

20 the probability values assigned to each braach are not normal- 
ized and can be greater than one. 

According to an inplementation mode of the above - 
described method, a weighting is assigned to each branch, the 
weighting of the branches of the last ramifications being equal 

25 to the product of the prc±>ability value assigned to this branch 
and of the number of combinations of this branch, the weighting 
of the other branches being equal to the sum of the weightings 
of the branches originating from the considered branch and being 
on the next ramification level. 

30 According to an inplementation mode of the above - 

described method, the probability value assigned to each branch 
can be nonnalized, the normalized probability value of a branch 
being obtained by dividing the probability value of this branch 
by the weighting assigned to the first branch of the tree. 
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According to an lirplementatlon mode of such a method, 
the selection of a coribination is performed by iitplementing a 
method generating combinations having high probability values. 

According to an iirplementation mode of such a method, 
5 the representation of the probability distribution of all the 
combinations is memorized and may be siabsequently refined by the 
creation of additional branches, or may be sinplified by the 
suppression of certain branches. 

According to an inplementation mode of such a method, 
10 the number of values likely to be taken by a parameter is 
artificially increased, the probability value of a cofribination 
of values of control parameters, among which at least a value of 
one of the parameters corresponds to an added value, is zero. 

The foregoing objects, features, and advantages of the 
15 present invention will be discussed in detail in the following 
nan- limiting description of specific embodiments in connection 
with the accoxipanying drawings, among which: 

Fig. 1 shows a car attempting to follow a truck; 

Fig. 2A illustrates a side view of the car of Fig. 1; 
20 Fig. 2B illustrates a top view of the car of Fig. 1; 

Fig. 3 illustrates a possible configuration of the car 
and of the truck as well as three futxire possible positions of 
the trucks- 
Figs. 4A to 4G illustrate steps of a method according 
25 to the present invention of construction of a representation of 
a probability distribution; 

Fig. 5 illustrates in the form of a ramified tree the 
steps of Figs. 4B to 4G; 

Figs. 6A to 6D illustrate two possible division modes 
30 of the set of couples (0^,0^) according to the method of the 
present invention; 

Fig. 7 illustrates an application of the method of the 
present invention to the failure diagnosis of an electric 
conponent assembly; and 
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Fig. 8 illustrates an application of the methcxi of the 
present invention to the recognition of figures. 

The method of the present invention applies to any 
system defined according to the criteria discussed hereafter. It 
5 can be decided of the value to be given to n specific parameters 
XSi to XSn of the system, knowing the values of n raeaisurement 
parameters XM^ to corresponding to a determined state of the 
system. Each of the parameters can take a finite number of 
values. The values of the specific parameters form a continuous 

10 sequence of integers. The measurement parameters may be symbolic 
variables, where the possible values can be yellow, blue, green, 
and red. A model of the joint probability distribution of the 
set of system parameters is known and the probability p(XSi,..., 
XSn,XM,...,XI^) for a given combination of all the parameters 

15 (XSi,...,XSn,XMi,...,XJ^) to be pertinent is calculated- The 
analytic fimctions and the databases used by the probability 
distribution model of the system are known. 

At any time, it is possible to make from the system 
model an inference consisting of defining the probability 

20 distribution of the coiribinations of values of all or part of the 
specific parameters, for exanple (XSi,XS2fXSn) # knowing the 
values of all or part of the measurement parameters, for 
exanple, (XMi,XM3). Probability p(XSi,XS2,XSn/XMi,XM3) for the 
selection of a combination (XSi,XS2,XSn) to be pertinent knowing 

25 the values of measurement parameters (XMi,XM3) is defined by: 

p(XS„XS2,XS„/XM„XM3) = i Xp(XS„...,XS„,XM„...,XM„) (3) 

XM2,XM4,...XM„ 

where Z is a normalization constant. 

After noting down the considered measvirement parame- 
ters, the present invention provides a method for constructing a 
30 representation of the probability distribution of the combina- 
tions of the values of the k selected specific parameters, 
obtained for the noted down values. A conbination of values of 
the k selected parameters will be called "a combination" hereaf- 
ter. The set of combinations may be represented by a set of 
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points defined in a k-dimensional space E. The probability 
distribution is then represented in a k+1 -dimensional space. 

The construction method of the present invention aims 
at dividing space E into several sets of points and at assigning 
5 an identical probability value to all the points of a same set 
to obtain a representation of the probability distribution of 
the combinations. 

Once the representation of the probability distribu- 
tion of the combinations has been obtained by the method of the 
10 present invention, the selection of a combination is performed 
according to one of several selection criteria. 

There can be a great variety of applications of the 
method of the present invention, as will appear from the reading 
of the mentioned exanples. 
15 In a first part, a method for constructing a represen- 

tation of the probability distribution of the combinations will 
be described. 

In a second psurt, different selection criteria will be 

described. 

20 In a third part, exanples of application of the method 

of the present invention will be described. 
1. CDNSTRUCTION METHOD 

1.1. General principle 

The method for constructing a representation of the 
25 probability distribution of the combinations consists of succes- 
sively selecting different combinations from among the set of 
possible combinations and of calculating their respective prob- 
ability values. After each selection of a combination, the set 
of points of space E containing the selected combination is 
30 divided into several sets of points. The set of points contain- 
ing the selected cOTibination takes the probability value of this 
cosnbination. The other sets of points keep the probability value 
that they had before the division. 
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Initially, space E is not divided and all the points 
of space E take probability value pi of the first selected 
conibination C^. 

The selection of a second conibination different 
5 from the first selected combination C^, caiises a division of 
space E into several sets of points, the set of points contain- 
ing the second selected conibination C2 taking probability value 
P2 of the second selected conibination €2, the other sets of 
points taking probability value pi of the first selected conibi- 
10 nation C^^. 

The selection of a third combination C3, different 
frcHu the first and second selected combinations and C2, 
causes the division of the set of "parent" points containing the 
third selected combination C3 into several sets of "child" 

15 points containing the third selected conibination C3 taking 
probability value P3 of the third selected conibination, the 
other sets of "child" points taking the probability value of the 
set of "pcorent" points, p^ or P2. 

The selection of a possible fourth conibination C4, 

20 different from those selected previously, would cause a new 
division of the set of "parent" points containing the fourth 
selected conibination C4 into several sets of "child" points, the 
set of "child" points containing the fourth selected combination 
C4 would then take probability value P4 of the foiirth selected 

25 combination C4, and the other sets of "child" points would take 
the probability value of the set of "parent" points pi, P2# or 
P3- 

This construction method can be repeated as many tiines 
as possible according to the available time. The representation 

30 of the probability distribution will be all the more accurate as 
the nxuiiber of selected combinations is large. Conversely to the 
creation of a representation according to the "gaussian mixture" 
method, the construction method of the present invention can be 
executed for a variable time, where the selection of the execu- 

35 tion time can be adapted to each system. 
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The successively-selected conibinatians can be obtained 
according to a pseudo-random method generating ccxnbinations 
xmiformly distributed over space E or according to an optimized 
method generating combinations having high probability values. 
5 According to an inplennentation mode of the method of 

the present invention, each set of "child" points, resulting 
from the division of a set of "parent" points, comprises an 
identical number of combinations. 

1.2. Construction tree 
10 The present invention provides keeping a trace of the 

construction of the probability distribution via a construction 
tree. 

The first branch of the construction tree represents 
space E and takes probability value Pi- The first branch 

15 ramifies into second branches each representing one of the sets 
of points resulting from the division of space E. Each second 
branch takes the probability value of the points of the second 
considered branch. 

The second branch associated with the set of "parent" 

20 points containing the third selected combination C3 ramifies 
into third branches each representing one of the sets of "child" 
points. Each third branch takes the probability value of the 
points of the third considered branch. 

Generally, ea:ch second branch is likely to ramify into 

25 several third branches according to the new selected combina- 
tions. Each third branch can ramify into several foxirth branches 
and so on. 

The construction tree is memorized along its construc- 
tion. The end branches of the tree give the final division of 

30 the set of combinations. The final representation of the pr<±>- 
ability distribution will be used hereafter to select one of the 
combinations, as will be described in the second portion. It may 
be provided to construct a representation of the probability 
distribution which is more or less accurate according to the 

35 available memory. 



The creation of a construction tree has severatL advan- 
tage, as will be specified hereafter, especially for the obtain- 
ing of nomnalized probability values and for the selection of a 
conibination by random drawing according to a random drawing 
5 method of the present invention. 

1.3, Illustration for the car/truck system 

1.3.1 Selection of a speed 
The construction of a representation of the probabil- 
ity distribution of the coaribinations according to the method of 
10 the present invention is illustrated hereafter for the previ- 
ously-described car/ truck system. 

The case where the car control system decides of the 
values to be given to control paLrameters {ol-^oX:' ^h) 
after the others is first considered. In the case of a speed 
15 selection, probability p(v) for the selection of a speed v at a 
time ti to be pertinent, knowing shift values c^i and Cyi noted 
down at time t^, can be calculated as follows: 

Fig. 4A shows a probability distribution 10 of the 
20 speed values obtained for shift values Cxq and Cyo noted down at 
time to* This distribution is that of which an approximation is 
desired to be obtained, siirply, rapidly, and minimizing the used 
calcxilation and memorization means. Speed v may take an integral 
value ranging between and including 0 and 15 km/hour. The speed 
25 is shown in abscissas, probability p(v) is shown in ordinates. 
In this exairple, probability distribution 10 of the speed values 
is a continuous function which is 0 when the speed is zero or 
greater than 14 and which exhibits two maximum values for speeds 
V equal to 4 km/h and 10 km/h. 
30 Figs. 4B to 46 altogether illustrate the construction 

of a representation of the probability distribution of Fig. 4A. 
The sets of speed values, or branches, are shown by a two-way 
corrow positioned imder the speed values belonging to the branch. 
A horizontal line cutting probability distribution 10, and 



placed above a two-way arrow, shows the probability vcd.ue asso- 
ciated with the speed values of the branch shown by the two-way 
currow. 

Fig. 4B illustrates a first step of the construction 
5 linked to the selection of a first speed vcdue equal to 4 
km/h of probability p^, where the set of speed values, forming a 
first branch B, takes probability value p^. 

Fig. 4C illustrates a second step of the construction 
linked to the selection of a second speed value V2 equal to 12 

10 km/h, of probability P2. The set of speed values is divided in 
two sets of speed values forming each of second branches and 
B2. Branch B^ gathers the smallest speed values ranging from 0 
to 7 km/h. Branch B2 gathers the highest speed values ranging 
from 8 to 15 km/h. The two branches B^ and B2 gather a same 

15 number of speed values. Branch B2 conteUJis the second selected 
speed value V2, it thus takes probability value P2. Branch B^ 
keeps probability value pi- 

According to a first aspect of the inplementation mode 
of the method of the present invention selected for this exam- 

20 pie, the ramification of a branch results in the creation of two 
branches conprising a same nuntber of speed values, one of the 
branches comprising the smallest speed values, the other branch 
coitprising the highest speed values. 

Figs. 4D and 4E illustrate two phases of a third step 

25 of the construction linked to the selection of a third speed 
value V3 equal to 6 km/h of probability value P3 . In a first 
phase, branch B^ containing the third selected speed value V3 
ramifies in two branches B^^^ and Bi^2 ^® appears in Fig. 4D. 
Branch B^^^ gathers the speed values ranging from 0 to 3 km/h, 

3 0 branch Bi^2 gathers the speed values ranging from 4 to 7 km/h. 
The first and third selected speed values v^ and V3 belong to 
the same branch Bx,2- this case, branches B^ and Bi^2 ^^^P 
probability value p]^. The ramification method will Ceunry on 
(Fig. 4E) until one of the branches only contains the third 

35 speed value V3. 
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According to a second aspect of the iitplementation 
mode of the method of the present invention selected for this 
exanple, the ramification of a "parent" branch containing the 
last selected speed value carries on until a "child" branch only 
5 contains the new selected speed value and no other previously- 
selected speed value. The intermediary "child" branches take the 
probability value of the "parent" branch. 

Fig. 4E illxistrates the second phase of the third 
step. Branch B^^2 containing the third selected speed value V3 

10 thus ramifies in two branches Bi^2,l ^1,2,2 ^respectively 

conprising speed values 4.5 and 6.7 km/h. Branch Bi 2,2 
contains the third selected speed value V3 and no other selected 
speed value. Probability value P3 is thus assigned to branch 
^1,2,2* Branch Bi^2,l keeps probability value p^ of branch Bi^2 

15 from which it originates. 

Fig. 4F illustrates a fourth step of the construction 
linked to the selection of a fourth speed value V4 equal to 10 
km/h of prtsbability value P4. Branch B2 containing the foiirth 
selected speed value V4 ramifies in two branches B2,i and 62^2* 

20 the branches respectively gathering speed values 8 to 11 and 12 
to 15 km/h. The fourth speed value V4 belongs to branch 62^1 and 
no other selected speed value belongs to this branch. Probabil- 
ity value P4 is then assigned to branch 62^1. Branch B2^2 keeps 
probability value P2- 

25 Fig. 4G illustrates a fifth step of the construction 

linked to the selection of a fifth speed value. V5 equal to 1 
km/h, of probability p^. Branch B±^i containing the fifth 
selected speed value V5 ramifies in two branches B^ ^ ^ and 
Bi,i,2' ^® branches respectively gathering speed values 0, 1 

30 and 2, 3 km/h. The fifth selected speed value V5 only belongs to 
branch B^^i^i and no other selected speed value belongs to this 
branch. Probability value P5 is then assigned to branch B^^^.,!' 
Branch B^ ^ 2 keeps probability value p^. It should be noted 
that at this last stage, a good approximation of probability 

35 distribution 10 of Fig. 4A has been obtained. 
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Fig. 5 shows the construction tree of the represen- 
tation of the probability distribution of the speed values 
obtained according to the five steps described in relation with 
Figs. 4A to 4G. Branch B of probability value pi ramifies in two 
5 branches Bi and B2 of respective probability values pi and P2. 
Branch B2 ramifies in two branches B2^i and B2^2 respective 
probabilities pi and P2. Branch B^ ramifies in two branches B^^i 
and Bi^2 of probability value px- Branch B^^i ramifies in two 
branches B^^i^i and Bi^i,2 respective probability values P5 

10 and pi. Branch Bi^2 ramifies in two branches Bi^2,l ^1,2,2 
of respective probability values pi and P3. The final division 
of the set of speed values is provided by the end branches of 
the construction tree. 

1.3.2. Selection of angles (ah,a^ r) 

15 The construction of a representation of the probabil- 

ity distribution of the combinations according to the method of 
the present invention is illustrated hereafter in the case where 
the car control system decides of the values to be given to two 
control parameters. 

20 Figs. 6A to 6D show the two-dimensional space E of all 

the coi^jles of values of horizontal and vertical inclination 
angles {0^,0^} . A couple of horizantatL and vertical inclination 
angles (ah/Ctv) will be called a couple hereafter. Horizontal 
inclination angle is shown in abscissas. Vertical inclination 

25 angle Oy. is shown in ordinates. 

Probability p(ah/CXv./cxi,Cyi) for the selection of a 
couple (ah/Oy.) to be pertinent knowing shift values Cxi and Cyi 
can be calculated according to the following formula: 

p(a„a,/c,j,Cyi) = -i: 5;p(v,a^,,a,,ah,c,i,Cyi) (4) 

30 where Z is a nonnalization constant and p(v,a^ot'^'^h'^'^y) 
calculated according to formula (1) . 

The implementation mode of the method of the present 
invention for this example uses the first and second above - 
described aspects. The branches resulting from a ramification 

35 are at the number of two and contain the same number of couples. 
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The ramification of a bz3nch carries on imtil the last selected 
cot^le is the only selected couple of one of the branches. 

Horizontal inclination angle aj^ can take six values 
between 0® and 5**, vertical inclination angle can take four 
5 values between 0<* and 3**. 

To sinplify the computer processing, the nxjihber of 
values likely to be taken by a parameter, if it is not a power 
of two, is increased to the immediately greater power of two. 
The probability of the couples having one of their parameters 
10 corresponding to an added value (not initially provided) is 
zero. 

In this exaitple, horizontal inclination angle can 
initially take six values. The number of values is thus artifi- 
cially brought to 8 (2^), and the possible values now cure from 
15 0® to 7®. No increase in the number of values is performed for 
vertical inclination angle for which 4 {2^) values are 
possible . 

Fig. 6A shows the set of couples {ai^,a^) forming first 
branch B. As previously, a first coi:?>le Ci (ahi=2, 0^1=3) and 

20 the probability value pi calculated for the first couple is 
assigned to all the couples of branch B. 

Fig. 6B illxistrates a possible ramification of branch 
B after selection of a second couple C2 (o^h2-'^' av2=^) 
probability P2. Branch B ramifies in two branches B^ and B2 

25 according to a vertical limit 12 passing between horizontal 
inclination angle values 3^ and 4®. Branch B^ gathers the 
couples having a horizontal inclination angle strictly smaller 
than 4 (2^) . Branch B2 gathers the couples having a horizontal 
inclination angle greater than 4 {2^) . 

30 AccorddLng to an aspect of the implementation mode of 

the method of the present invention for this exairple, the 
ramification of a "parent" branch results in the creation of two 
branches according to a vertical limit. In the case where it is 
impossible to define a vertical limit, that is, when the couples 

35 of the "parent" branch all have the same horizontal inclination 



value tthf the division is perfoimed according to a horizontal 
limit passing between two vertical inclination values a^. 

Figs. 6C and 6D illustrate another possible ramifica- 
tion of branch the second selected couple C21 (ay2i=lf 
5 ah2'=l) being different frcm C2. Pig. 6C illustrates a first 
division of branch B according to the same vertical limit 12 as 
that previously defined. The selected first coi?>le Ci and second 
couple C21 both belong to branch B^. Branches B^ and B2 take 
probability value Pi of first couple and a new ramification 

10 of branch B^ containing couple C21 is performed until the two 
selected couples and C21 are in different branches. 

Fig. 6D illustrates the ramification of branch B^ 
according to a horizontal limit 13 passing between values 1° and 
2^ of vertical inclination angle a^. Branch B^^i gathers the 

15 couples having a vertical inclination angle value greater than 
or equal to 2. Branch Bi^2 gathers the couples having a vertical 
inclination value strictly smaller than 2. Branch Bi^2 take 
probability value P2 of second couple €2* and branch B2,2 takes 
probability value p^. 

20 According to another aspect of the implementation mode 

of the method of the present invention for this exaitple, the 
successive ramifications of a "parent" branch are successively 
performed according to a vertical limit and a horizontal limit 
xantil the new selected cartibination is the only one of the 

25 selected cOTibinations to belong to a given "child" branch. 
1.4. Normalized probability values 

The probability values, for exaitple p(XSi, XS2/ XS^) , 
obtciined from formula (3) are in fact defined to within a 
constant, normalization constant Z. This constant can be calcu- 
30 lated only when the probability values of all canbinations (XS^, 
XS2, XSn) are known and have been calculated without taking Z 
into account (Z taken to be equal to 1) . Normalization constant 
Z is then equal to the sum of all the probability values. 

In practice, only very few probability values are 
35 calculated on construction of the representation of the 
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probability distributicm. The prcibability values calculated for 
the selected ccnibinations are not normalized. 

The present invention provides defining an interme- 
diary normalization constant Zi which is an estimate of normali- 
5 zation constant Z. Normalization constant Zi is calculated 
during the construction of the probability distribution 
representation. It is equal to the sum of the probability values 
of all corabinatians, the probability value of a carbination 
being that assigned to the branch containing the considered 

10 combination. 

To easily calculate intermediary nortnalization 
constant Zi, the present invention provides assigning a weight 
to each branch. The weighting of the branches of the last rami- 
fication is equal to the product of the probability value 

15 associated to the considered branch and of the number of combi- 
nations of this branch. The weighting of the other branches is 
equal to the svrm of the weightings of the branches originating 
from the considered branch and located on the next ramification 
level. The weighting of the first branch is then equal to inter- 

20 mediary normalization constant Zi. 

The weightings are updated on ramification of one of 
the end branches of the construction tree. The weightings of the 
branches located between the first branch and the ramifying end 
branch must be recalculated. 

25 Thus, for the example of construction of the speed 

probability distribution representation shown in Figs. 4B to 4G, 
the weighting of branch at the end of the first step, Fig. 
4B, is 16*pi. At the end of the second step. Fig. 4C, the 
weightings of the new branches B^ and B2 respectively are 8*pi 

3 0 and 8*P2i and the weighting of branch B is updated and now is 
8*pi+8*P2. At the end of the third step. Fig. 4D, the weightings 
of branches Bi^i and Bi^2 both are 4*pi, and the weighting of 
branches B^ and B remains unchanged (respectively 4*pi+4*pi=8*pi 
and (4*pi+4*pi)+8*pi=16*pi) since the probability values 

35 assigned to the different speeds have not changed, only the 
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division of the speed values having cdianged. At the end of the 
fourth step, Pig. 4E, the weightings of branches Bi^2,l 
^1,2,2 respectively are 2*pi and 2*P3, the weighting of branch 
^1,2 is 2*pi+2*P3, the weighting of branch is 

5 4*pi+(2*pi+2*P3) and the weighting of branch B is 
(4*pi+ (2*Pi+2*P3) ) +8*P2 . 

An advantage of the method of the present invention is 
that it enables rapidly knowing the nomnalized probability value 
of the speed since it is not necessary to calcxilate all the 
10 probability values. 

1.5. Advantages 

The method of the present invention enables represent- 
ing a great variety of probability distributions, conversely to 
the so-called "gaussian mixtvire" method. 

15 The method of the present invention enables obtaining 

a more or less detailed representation according to the avail- 
able memory and to the allowed calculation time. 

Further, the iitplementation of an optimized method for 
the combination selection enables obtaining within a very short 

20 time a representation taking vp little memory and exhibiting a 
sufficient accuracy for the areas of space E where the combina- 
tions exhibit high probability values. The method of the present 
invention thxis is "multi-resolution" in that the division of 
space E can be very fine for certain portions and rough for 

25 others. 

This "multi-resolution" feature of the method of the 
present invention enables, conversely to existing representa- 
tions, constructing within a relatively short time and by using 
memories of reasonable size, probability distribution represen- 
30 tations of a large number of control parameters or parameters 
that can take a large number of values. 
2. SELECTION 

Kiiown modes of selection of one of the ccxiibinations 
based on a representation of the probability distribution of the 
35 cambinations consist of selecting one of the combinations 
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exhibiting the naximum probability or of selecting a combination 
by random drawing, the principle of which will be recalled 
hereafter. Other selection criteria such as the selection of the 
one of the combinations having a predetermined probability value 
5 or the selection of one of the cosiibinatians having a probability 
value ranging between two given probability values may however 
be defined. 

2.1. Selection of a combination of maximum probability 
According to an inplementation irKxie of the method of 

10 the present invention, a memory register is used to store an 
indication of the branch (es) exhibiting the maximum probability 
value- The register initially n^morizes the first branch of the 
construction tree, after which it is updated during the 
construction of the probability distribution representation. At 

15 each ramification of a branch, it is checked whether the 
probability value of the last selected combination is greater 
than the probability value of the branch memorized by the regis- 
ter, and if such is the case, the register is updated and 
memorizes the new branch containing the last selected conibi- 

20 nation. 

In the case where the branch of maximum probability 
ramifies and gives one or several "child" branches of same 
probability, the register is updated and memorizes all the 
"child" branches. 

25 The selection of a combination then consists of 

identifying the branch memorized by the register containing the 
greater number of combinations and of then selecting one of the 
combinations of this branch. 

2.2. Selection by random drawing 

30 A random drawing consists of selecting one of the 

possible combinations in a way such that a combination exhibit- 
ing a high probability stands in a good chance of being selected 
and a combination exhibiting a low probability stands in a poor 
chance of being selected. After a great number of random draw- 

35 ings, the probability distribution of the "drawn" combinations 



is identical to the probability distribution of the initial 
conibinations on which the randcm drawing method is based. 

According to an inplementation mode of the method of 
the present invention, the random drawing of a ccnbination is 
5 performed according to a recursive selection method using the 
construction tree of the representation of the probability 
distribution . 

Starting from the first branch, a second branch, then 
a thiixi branch are selected, and so on \mtil an end branch of 

10 the construction tree is selected. For this purpose, a second 
branch is selected by performing a random drawing from among the 
second branches. The second branches exhibiting the highest 
probabilities stand in the best chance of being selected, and 
conversely. One of the third branches originatdLng from the 

15 selected second branch is similarly selected by performing a 
random drawing between these third branches, and so on. 

In the case where each branch ramifies in two 
branches, the method of the present invention conprises several 
steps described hereafter. 

20 In a first step of the method, a nxmiber p ranging 

between 0 and 1 included is randomly selected. 

In a second step, sum S of the probability values 
assigned to the 2 "child" branches resulting from the ramifica- 
tion of the first branch is calculated- Then, for each "child" 

25 branch, a new probability value equal to the ratio between the 
probability value assigned to the considered branch and the 
calculated sum S is calculated. 

In a third step, two contiguous probability intervals 
are defined between 0 and 1, the first interval being associated 

30 with a first child branch, the second interval being associated 
with the second child branch. The first interval ranges from 0 
to and including the probability value of the first child branch 
and the second interval ranges from this probability value to 1. 




In a fourth step, it is identified in which interval 
nvmiber p is located and the "child" branch eissociated with this 
interval is selected. 

In the case where the selected "child" branch ramifies 
5 in other branches, the recursive method is resumed at the first 
step, the first branch being replaced with the selected "child" 
branch. The first selected number p may be used again. 

In the case where the selected "child" branch does not 
ramify into other branches, the recursive method stops and one 
10 of the combinations of the selected child branch is selected. 

The above-mentioned recursive selection method is used 
in the exaitple of the car/truck system to select a speed or a 
couple of inclination angles {a^,ai^) . 

An advantage of the method of the present invention is 
15 that the random drawing method is sinple and easy to inplement. 

2.3. Tree memorization 

Once the decision has been taken, the construction 
tree can be erased from the memory. It may however be provided, 
for systems such as the car and the truck, to keep in a "cache" 

20 memory the construction trees successively obtained after 
different samplings of the measurement parameters. The most 
often used construction trees may be memorized. 

In the case where a probability distribution is often 
used, its representation can be refined by carrying on the 

25 ramification of the construction tree of this representation. In 
the car/ truck exairple, the tree ramification may be carried on 
for the time allowed between the sanpling of the values of c^ 
and Cy and the time when the values of a^ot' ^h nnist 

be selected. 

30 Similarly, in the case where a probability distribu- 

tion has been much used at a given time, and less afterwcuxls, 
the accuracy of the probability distribution can be decreased by 
suppressing more or less end branches of the construction tree. 

In the case where it is chosen to calculate a 

35 weighting for each branch to know the intermediary normalization 
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constant such as defined hereabove, the weighting values of the 
branches located between the first branch and the suppressed 
branch (es) will be updated. 

An advantage of the method of the present invention is 
5 that it enables refining or siirplifying a probability distribu- 
tion representation. This enables implementing strategies of 
memorization of different probability distributions to globally 
iitprove the successive combinations selections. 
3. EXAMPLES OF APPLICATION 

10 3.1. Failure diagnosis 

Fig. 7 shows an electric device that conprises several 
cortponents between an input I and an output O, each component 
being capable of conducting part of input current K when operat- 
ing, with no current flowing when the coirponent is defective. 

15 A conponent A is placed between input I and a first 

intermediary point 100. Output current KA of conponent A is 
equal to 100% of current K when the conponent operates (and 
equal to 0% of current K when the conponent is defective) . 

Conponents B and C are placed in parallel between 

20 first intermediary point 100 and a second intermediary point 
101. Output current KB of conponent B is at most equal to 40% of 
current K when conponent B operates. 

Conponent C is formed of two conponents Ci and C2 in 
parallel. Each conponent Ci and C2 can conduct up to 30% of 

25 current K. Output current KC of conponent C may thus be equal to 
0%, 30%, or 60% of current according to whether the two 
ccxrponents are defective or the two conponents are operative. 

A conponent D is placed between point 101 and output 
O. Conponent D is formed of eight conponents to Ds in paral- 

30 lei. Each COTponent to Ds can conduct up to 15% of current K 
when it operates. Further, it is necessary that at least six of 
conponents Di to Ds operate for conponent D to operate . 

Current KBC entering conponent D is equal to the sum 
of c\arrent KB and of current KC. Current KBC may be equal to 0%, 

35 30% (only C^ or C2 is operative), 40% (only B operates), 60% (C^ 
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and C2 are operative) , 70% (Ci or C2 and B are operative) or 
100% (Ci, C2, and B operate) of current K. 

Output current KD of conponent D may be equal to 0%, 
30%, 40%, 60%, 70% (KBC is respectively equal to 0%, 30%, 40%, 
5 60%, and 70% of current K and at least six of conponents Di to 
Ds cure operative) , 90% (KBC is equal to 100% of current K and 6 
conponents to D3 are operative) or 100% (KBC is equal to 100% 
of current K and at leaist 7 conponents to D3 are operative) 
or current K. 

10 Each conponent A, B, C^, C2 and to Dg can be 

considered as being a parameter of the device that can take 
value 0 when the ccxtponent is defective, and 1 when the conpo- 
nent is operative. Output currents KA, KB, KC, KBC, and KD are 
also considered as being parameters of the device that can take 

15 more or less values. Thus, KA and KB can take value 0 or 1 
according to whether the current is respectively 0% or 100% of 
current K. KC can take values 0, 1, or 2 according to whether 
the current respectively is 0%, 30%, or 60% or current K. KBC 
can take values 0, 1, 2, 3, 4, or 5 according to whether the 

20 current respectively amoimts to 0, 30, 40, 60, 70, or 100% of 
current K. KD can take values 0 to 7 according to whether the 
current respectively is 0, 30, 40, 60, 70, 90, or 100% of 
current K. 

Probability p(A,B,Ci,C2/Di,...,D8,KA,KB,KC,KBC,KD) for a 
25 given combination of all parameters to be pertinent is defined 
as follows: 

p(A,B,Ci,C2,Di,..„D8,KA,KB,KC,KBC,KD) = 
p(A)p(B)p(Ci)p(C2)p(Di)p(D2)p(D3)p(D4)p(D5)p(D6)p(D7)p(D8) 
p(KA)p(KB)p(KC)p(KBC)p(KD)p(KA/A)p(KB/B,KA)p(KC/C,KA) 
30 p(KBC/KB,KC)p(KD/KBC,D) (4) 
where 

p(A), p(B), p(Ci), p(C2), p(Di) to p(D8), p(KA), p(KB), p(KC) , 
p(KBC), and p(KD) respectively are the probabilities for A, B, 
Ci, C2/ KA, KB, KC, KBC, and KE to have a given value. 
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P(KA/A) is the probability for KA to have a given value knowing 
the value of A, 

p(KB/B,KA) is the probability for KB to have a given value 
knowing the value of B and of KA, 
5 p(KC/C,KA) is the probability for KC to have a given value 
knowing the value of C and of KA, 

p(KBC/KB,KC) is the probability for KBC to have a given value 
knowing the value of KB and of KC, 

P(KD/KBC,D) is the prc±ability for KD to have a given value 
10 knowing the value of KBC and of D, 

Siitple probabilities p(A) to p(KD) can be defined frcxn 
databases noting down the failure cases having appeared in tests 
performed by the firm manufacturing the device. Conditional 
probabilities p(KA/A) to p(KD/KBC,D) can easily be defined based 
15 on the device description. For exaitple: 

p(KA^O%/A;=0)=l 

p(KA=100%/A=0)=0 

p(KA=O%/A=l)=0 

p(KA=100%/A=l)=l 

20 In a failure diagnosis of the device, the ciunrent at a 

point of the device (KA, KB, KC, KBC, or KD) can be measxired 
and/or the operating state of one or several corrponents of the 
device (A, B, C^, C2, to Ds) can be considered. Once one or 
several currents have been saitpled and/or one or several coitpo- 

25 nents have been analyzed, it can be determined which coirponents 
are likely to fail (diagnosis 1) or what current is likely to 
flow at a point of the device (diagnosis 2) . 

According to the performed diagnosis, parameters A, B, 
Ci, C2, T>i to Ds, KA, KB, KC, KBC, and KC can be either specific 

30 parameters, the value of which is desired to be determined, or 
measu3rement pcurameters, since the current or the operating state 
is measured, or non- retained parameters since their value is not 
desired to be determined and their value is not measured. 

In diagnosis 1, the measurement parameters are one or 

35 several currents, for exanple, KBC, and the specific parameters 
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are one or several ccnponents, In this exanple, conponents A, B, 
Ci, and C2 located i5)streain of current KBC. After a saitpling of 
the value of current KBC, a representation of the probability 
distribution of the combinations of values of parameters A, B, 
5 C^, and C2 is constructed according to the method of the present 
invention, knowing the value of measiireitient parameter KBC. 
Probability p(A,B,Ci,C2/KBC) for a corribination (A,B,Ci,C2) to 
represent the state of the device, knowing the value of KBC, can 
be calculated as previously by summing up all the prc±>abilities 

10 p(A,B,Ci,C2,Di,_,D8,KA,KB,KC,KBC,KE) for all the parameter 
values not retained in this inference (Di to Dg, KA, KB, KC, 
KD) . It will be first desired to identify a first ccmbinatian, 
that exhibiting the maximimi probability, to perform a first 
repair. If this repair appears to be insufficient or unfounded, 

15 a second conibination exhibiting the highest probability after 
the first combination will be identified, and so on. 

In diagnosis 2, the control parameter will be the 
CTirrent which is desired to be known, for exanple, KBC, where 
the measurement parameters may be one or severed, of the other 

20 parameters. After a sanpling of the vcQue of the meeusurement 
parameters, a representation of the probability distribution of 
the possible values of current KBC will be constructed according 
to the method of the present invention. The value of current KBC 
exhibiting the maximum probability is then selected. 

25 3.2. Figure recognition 

Pig. 8 shows an image 200 of a figure written by hand 
which is desired to be identified. The image is broken down into 
64 squares or pixels of identical sizes. For each pixel, a grey 
level representing the surface area taken up by the lines of the 

3 0 figure in the considered pixel is measured on a scale from 0 to 
16. 

The system carprises a single psurameter to be deter- 
mined (or specific parameter) CHIFFRE that can take values 0 to 
9 and 64 meaisurement parameters Pix [i] , i being an integer 
35 ranging between 1 and 64, that can each take values 0 to 15. 
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Probability p(CHIFFRE,Pix[l] ,«,Pix[64] ) for a given 
combination of values of all the parameters to be possible, is 
defined as follows: 

p (CHIFFRE , Pix [1] , ^, Pix [64] ) = 
5 p (CHIFFRE) *p (Pix [1] /CHIFFRE) *^*p (Pix [64] /CHIFFRE) (5) 
wliere p (CHIFFRE) is tlie probability to liave a given figttre, and 
p (Pix [i] /CHIFFRE) is the probability to have a given grey level 
for pixel Pix [i] , knowing the figure . 

It is generally considered that figures 0 to 9 are 
10 equiprobable and thus that p (CHIFFRE) is equal to l/lO. Condi- 
tional probabilities p (Pix [i] /CHIFFRE) are defined at the end of 
a training phase consisting of measuring the grey levels of each 
pixel for different models of handwritten figures. Each prob- 
ability p (Pix [i] /CHIFFRE) can be defined by a histogram (with 16 
15 columns) normalized according to the Laplace law. 

The recognition of a figure conprises a first pliase of 
measurement of tlie grey level of each pixel. In a second pliase, 
a representation of the probability distribution of the figures 
is constructed l>ased on formula (5) and on tlie method of the 
20 present invention, Imowing the grey levels of each pixel. The 
selection of a figure consists of identifying the figure exhib- 
iting the maximum probability. 

3.3, Evaluation of a transport cost 

A shipping COTpany dispatches by cargo goods from a 
25 European port to another port. The dispatched goods can be of a 
great variety: food products, medicine, electric appliances, or 
clothes. Different sorts of containers are used for the storage 
of the goods in the cargo and in the port of arrival. The 
containers may be refrigerating and of different sizes. 
30 The shipping coitpany desires to rapidly determine, in 

a telephone conversation, for exairple, transport expense esti- 
mates Imowing the port of departxire (PortDep) and the port of 
arrival (PortArr) , the type of transported goods (mr) , the 
container used (Cant) , the client (Cl) , and the month (M) in 
35 which the transport will occur. These parameters, previously 
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input on deterrnination of the estimate, form the measiarement 
parameters of the transport system. 

Further, the ccnpany has in its possession a whole set 
of information especially relative to the contcdner prepcoration 
5 time in the European port of departure (TdP) , the shipping time 
(TdlM) (outward or back, the cantsmiers are full on the outward 
trip and enpty on the way back) , the waiting time of the 
container in the port of arrival (TdA) , the time of container 
\mloading at the client ' s (TdDC) , the container reconditioning 
10 time when back in Europe (TdRE) , the times being counted in 
days. 

Similarly, information is available relating to the 
daily cost of the renting of a container (CdLJ) , the daily cost 
of immobilization of a container for its reconditioning (CdlJ) , 

15 the shipping cost {Cd™) , the repair cost of a container (CdR) . 

The total transport cost (CT) is the sum of the total 
cost of the container renting (CdLT = CdLJ* (TdP + 2*TdTM + TdA + 
TdDC + TdRE) ) , of the totcLL container immobilizatian cost 
(CdIT), of the shipping cost (CdTM = CdIT*TdRE) , of the cost of 

20 the load balancing in the cargo (CdE) , and of the cost of a 
container repair (CdR) . 

The mcxJel of the joint probability distribution of all 
the previo\isly-de fined parameters is constructed based on the 
independent probability distributions defined for each time 

25 pcurameter (TdP, TdTM, TdA, TdDC, and TdRE) , and for each cost 
parameter (CdLT, CdIT, CdTM, CdE, CdR, CT) . The probability 
values are obtained from a set of data acquired along the 
coitpany's lifetime. For exanple, the probability distributions 
of the container preparation time in the port of departure (TdP) 

30 knowing the nature of the good (Mar) are a gaussian family. The 
probability distributions of the shipping costs (CdTM) knowing 
the type of container (Cent) , the port of departxire (PortDep) , 
the port of arrival (PortArr) , and the transported goods (Mar) 
are Dirac functions. 
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Generally, the total transport cost (CT) is desired to 
be defined without detailing the intermediary costs. In this 
case, there is a single parameter (or specific parameter) to 
determine : the total cost (CT) . To rapidly provide a total 
5 cost, a representation of the prdbability distribution of the 
possible total costs CT is constructed according to the method 
of the present invention, knowing the measurement parameters 
(PortDep, PortArr, Mar, Cont, CI, and M) . The vendor first wants 
to know what the maximum cost is, for exairple, 2,000 eviros. He 

10 then establishes a first estimate by possibly taking a 10% 
margin with respect to the maximum cost, then offers 2,200 
euros. In the case where the client does not accept this price, 
the vendor can then estimate the average cost or the cost range 
containing for exanple 90% of the possible cost values. The 

15 average cost can easily be calculated by dividing intermedicury 
normalization constant Zi by the nioniber of possible total costs. 
The vendor will then establish a second estimate by taking a 
margin that may be smaller with respect to the average cost or 
with respect to one of the costs of the cost range. 

20 The established estimate may also detail all the 

costs, in which case the measurement control parameters of the 
system are (CdLT, CdIT, CdTM, CdE, CdR) • The total cost is then 
calculated based on the retained cost combination. 

This exanple of application shows that based on a 

25 representation, the combinatian of maximum probability can 
easily be determined, the average value and the standard devia- 
tion of the values of a parameter can be calculated, and one of 
the combinations exhibiting a given probability can be selected. 
The method of the present invention thus enables easily iirple- 

3 0 menting several selection criteria. 

Of course, the present invention is likely to have 
various alterations, modifications, and improvements which will 
readily occur to those skilled in the art. In particular, it 
will be within the abilities of those skilled in the art to 

35 define the branch ramification method most adapted to the 
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Studied system. Similarly, it will be within the abilities of 
those skilled in the art to define new criteria of the selection 
of a ccmnbination based on the tree-shaped representation of the 
probability distribution. 



